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Surface term for the capillary condensation transitions in a slit geometry 
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It is shown that a bare simple fluid model (SFM) proposed some years ago for studying adsorption 
between two semi-infinite solid walls can be improved by modifying the surface term in the grand 
potential for the film phase. Such a correction substantially improves the agreement between the 
predictions for phase transitions provided by that SFM and results obtained from calculations carried 
out for ^He with the density-functional method at zero temperature. The corrective term depends 
on the strength of the adsorption potential and observables of bulk helium. 

PACS numbers: 61.20.-p, 64.70.-p, 67.70.+n 
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A simple fluid model (SFM) has been proposed by 
Gatica et al. Q for exploring the behavior of adsorp- 
tion between two parallel walls separated by a distance 
L. In that paper, transitions between empty (E), film 
(F), and capillary condensation (CC) phases were deter- 
mined. Subsequently, Calbi et al. have checked such a 
SFM by comparing its predictions with results provided 
by the nonlocal density-functional (NLDF) formulated 
by the Orsay-Paris (OP) collaboration 3] for superfluid 
helium. By looking at Figs. 15 and 16 in Ref. 2i one may 
realize that the E ^ F phase transitions given by SFM 
for ''He/Li and ^He/Au systems do not reproduce quite 
well results obtained from the OP-NLDF theory. In this 
report we suggest a correction which substantially im- 
proves the agreement. 

For the descriptions provided in Refs. 1 and 2 it is 
assumed that ''^He atoms interact with the walls via the 
net potential 

U,^it{z)^V{z) + V{L-z), (1) 

with V{z) being the standard (9,3) adsorption potential 

V{z) = [ACll{27D^)]z-^ - , (2) 

where D is the well depth and C3 is the strength of the 
asymptotic van der Waals interaction |^]. This potential 
V{z) exhibits a minimum at 



[2C3/(3i^)]i/3 



(3) 



A sample of parameters D and C3 corresponding to the 
interaction of ^He atoms with several substrates given in 
Refs. 1^ and 1^ is listed in Table U The reader may find a 
plot of D vs C3 in Fig. 2 of Ref. 5. It should be mentioned 
that for alkali metals the adsorption potentials of Ref. ^ 
are more attractive than the previous ones of Ref.H 



TABLE I: Values of the weU depth D, the van der Waals 
strength C3, the position of the potential's minimum Zm, the 
reduced well depth D* and other relevant quantities. 
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jja 






D' 


a* 


/?PH 




/3mho 




[K] 


[kA3] 


[A] 












Na" 


12.53 


1197.4 


3.99 


8.03 


0.2177 




2.65 


2.18 


Li' 


17.87 


1422.5 


3.76 


10.78 


0.2054 


2.80 


2.81 


2.53 


Ha 


28.00 


360.0 


2.05 


9.20 


0.2487 


2.23 


2.32 


2.34 


Mg' 


35.63 


1850.2 


3.26 


18.64 


0.1856 


3.12 


3.11 


3.32 


Ne 


60.00 


163.0 


1.22 


11.74 


0.2664 


2.78 


2.16 


2.64 


Al 


60.30 


2340.0 


2.96 


28.63 


0.1708 


3.77 


3.38 


4.12 


Cu 


63.80 


2610.0 


3.01 


30.83 


0.1669 


3.90 


3.46 


4.28 


Ag 


63.80 


2890.0 


3.11 


31.90 


0.1641 


3.95 


3.52 


4.35 


NaCl 


69.60 


1230.0 


2.28 


25.42 


0.1879 


3.73 


3.07 


3.88 


Ar 


78.90 


870.0 


1.94 


24.63 


0.1970 


3.83 


2.93 


3.82 


Au 


92.80 


3180.0 


2.84 


42.27 


0.1566 


4.59 


3.69 


5.01 


LiF 


94.00 


1080.0 


1.97 


29.75 


0.1872 


4.26 


3.08 


4.20 


Gr 


190.00 


2130.0 


1.96 


59.64 


0.1577 


6.63 


3.66 


5.95 



"Potential parameters D and C3 are taken from Ref. 
except where noted. 'Potential parameters from Ref. 



According to SFM, the E ^ F transition for films of 
thickness £ adsorbed on both walls of the slit is governed 
by the evolution of the reduced (ie., in units of the sur- 
face tension cr/^ = 0.272 KA~^ [3) grand potential per 
unit area given by Eq. (8) of Ref. [l| 



A~D* [g{l + x] + g{L* - 1} 
-g{L*-l-a;}]+2a;A. 



(4) 



Here the first term is the surface energy, the second is due 
to the helium- walls interaction, and the last one measures 
the departure from the saturated vapor pressure (SVP) 
condition. The reduced well depth and gap are 



D* 



2pnD z„i/(7iv 



(5) 



*Also at Carrera del Investigador Cienti'fico of the Consejo Nacional 
de Investigaciones Cienti'ficas y Tecnicas, Av. Rivadavia 1917, RA— 
1033 Buenos Aires, Argentina 



A = (^0 - m) Po Zni/ai^ , (6) 

where pa = 0.021836 is the saturation equilibrium 
density and fig = —7.15 K the corresponding chemical 
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FIG. 1: Reduced phase diagram for *He confined by two pla- 
nar walls of Li (D* — 10.78, Zm = 3.76 A). Solid curves 
are SFM predictions and full circles OP-NLDF results, both 
of them calculated in the present work (see text). The 
dashed curve shows how much the correction introduced in 
the present work improves the prediction of the SFM for the 
E — > F phase transition. 

potential Pj. Furthermore, 

r , . r-U^AZrnZ*). 11 3 1 1 1 

^^^^ = A [ D ]^ 16-4^+167' 

here all the reduced distances are given in terms of Zm 
z* = zjzm, X = e/zm , and L* = L/Zm , (8) 

where x is the film thickness that minimizes ftp. 

Figure n shows that according to present calculations 
for the case of ^He confined into slits of Li the reduced 
grand potential flp given by the SFM becomes negative 
for i* > 9 at AsFM = 0.73, while the OP-NLDF predicts 
that the E — > F phase transition occurs at Adf = 0.85. 
Similar effects are also obtained for that substrates listed 
m Table |I| which are strong enough to produce a stable 
monolayer film. In order to illustrate this feature results 
for ■^Re/Mg and "^He/Au are plotted in Figs. Handel 
Let us recall that a slit of Na is too weak to produce a 
stable monolayer structure. 

In Figs. I1I3I besides the E F phase transition there 
are also indicated the E CC and F — > CC ones. These 
phase transitions are determined from comparisons with 
the corresponding grand potential per unit area obtained 
from Eq. (5) of Ref . il|, i.e., 

f^cc = — = 2-Z?*5{L*-l} + (i*-2)A. (9) 

In this work we shall demonstrate that the prediction 
given by the bare SFM is improved by introducing a cor- 
rection to the surface term in the grand potential 
given by Eq. An important shortcoming is that this 
term does not vanish in the limit of a zero-thickness film 
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FIG. 2: Same as Fig. but for *He confined by Mg {D* = 
18.64, Zm = 3.26 A). 
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FIG. 3: Same as Fig. d but for *He confined by An (D* = 
42.27, Zm. = 2.84 A). 

(i.e. for X — £/zm ^0). In order to eliminate this fail- 
ure one may follow the idea adopted by Cheng et al.^ 
in writing their Eq. (2.4). So that, we shall assume that 
the surface contribution grows exponentially from zero at 
£ — to the bare value 4 over a characteristic length ( 

17^(surf) = 4[l-exp(-£/C)] =4[l-exp(-/3a;)], (10) 

with 

/3 = z„JC = 1/C • (11) 

The parameter f3 was determined for all the systems 
quoted in Table HI by requiring that OP-NLDF results be 
reproduced by SFM. This means, we impose that the E — s- 
F transition must occur at Adf for both models. These 
phenomenological values denoted /3pH are also included 
in Table U The changes of the energetics due to the cut- 
off factor introduced in Eq. H10() only become important 
for rather thin films. 

It should be also noted that, if a similar correction 
is included in the grand potential for the CC phase the 
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corrected curves for the E — > CC and F — > CC transitions 
do not differ much from that displayed in Figs. 

The following lines are devoted to propose a possible 
origin for the cut-off factor /3 = Zm/C- As a first step, 
it is plausible to assume that the characteristic length Q 
be related to the width of the ground-state (g.s.) wave- 
function (submonolayer solution) of a simple harmonic- 
oscillator approximation to the adsorption potential in 
the neighborhood of each wall 



Kpp(2) = -D + Vno{z) ^-D+^k{z~ z„,f 



with the force constant 



k 



18 



3D 
2Ch 



5/3 



27- 



D 



(12) 



(13) 



see comment quoted as citation 44 in Ref. 
sionless force constant becomes 

K ^kzl^/D = 27. 

So, the reduced version of Va-pp{z) reads 



The dimen- 



(14) 



app 



^ ^app(gn^*) ^ _^ ^ ^ _ 



(15) 



and the reduced hamiltonian may be conveniently writ- 
ten as 



-"ho 



H* 



1 = -- 



2TOeff dz*"^ 

with an effective helium mass 

mcff ^mD z^ 
The g.s. energy of H* is 



f (16) 



(17) 



-^app - -1 



1 h 



2D 



(18) 

where = /rrics is a dimensionless kinetic energy fac- 
tor. Let us mention that as quoted in Ref. ,9^ for graphite 
(Gr) the semi-empirical He potential of Carlos and Cole 
[1(TI |. based on He scattering data, has k = 26. 
The probability density of the g.s. solution of Eq. 

is 



1 



/27ra* 



with 



m k 



27 



1/4 



(19) 



(20) 



FigureQlshows examples of density profiles yielded by the 
DF theory for the lowest stable coverages. Let us recall 
that in the literature the surface thickness is defined as 
the distance over which the density falls off from the 90% 
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FIG. 4: Comparison of monolayer *He systems adsorbed on 
the left wall of broad slits of Mg, Al, and Au. These films 
were determined at the E — > F phase transition, i.e., they 
correspond to the lowest stable coverage in each case. 



to 10% of the saturation value |T3|. In the case of the 
g.s. wavefunction given by Eq. ((T^ these quota lead to 



7^(C 



nCf) 



V'o(i + C) 

I V'o(l) P 



^^0(1 + CP 



0.90. 



(21) 



From here one gets 



/-* 



4 In 3 



c; + c 

= 1.69(7* 



(a*) 



= 0.10. (22) 



2x/21n3 



\/lnlO+ v/ln(10/9) 
\/3ct* . (23) 



It is interesting to notice that by assuming C,* 
the integral 



AC) 



1 



/2Tia* 



i+C* 



dz* e 2 



0.90, 



(24) 

amounts about 90% of the total norm. 

Turning to the cut-off parameter /3, the expression de- 
rived from a pure harmonic oscillator analysis becomes 



1 



1 

71 



27 



1/4 



1/4 



(25) 



The estimations provided by this formula are listed in 
Table From a glance at that table, one can realize 
that there is a very good agreement for the most weakly 
bounded monolayer systems, i.e., that confined by walls 
of H2 , Li and Mg. It is worthwhile to mention that the 
values of /? for ''He/Li and "^He/Mg are in the parameter 
region utilized by Cheng et al. for showing how SFM 
results approach those given by a full NLDF theory in the 
case of adsorption on semi-infinite surfaces with planar 
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the total effect be proportional to 



FIG. 5: Parameter f3pu obtained from equating SFM and DP 
results for ''He films as a function of \/ D* . Pull circles are 
metallic substrates, while open triangles stand for other kind 
of materials. The star is the /3mho value for a slit of Na. The 
solid straight line is given by Eq. CT . 



geometry. In particular, /3 = 3 is just the value chosen in 
that paper to illustrate the behavior. 

The values of /3ho fall short if the helium atoms are 
exposed to stronger attractions. Hence, one could ar- 
gue that a correction proportional to D should be intro- 
duced if the harmonic oscillator approach for the poten- 
tial V{z) of Eq. (|2Jl yields a g.s. oscillator energy much 
smaller than the well depth, i.e., when in Eq. H18|) holds 

^1. In searching for such a dependence, it was 



1 _fi 

2 D 

realized that the values of /3pH become aligned when plot- 
ted as a function of -\/ D* . This behavior may be observed 
in Fig. |5l where Na indicates the threshold obtained for 
the formation of stable monolayer films. Similar results 
are obtained by using the more recent NLDF of Ref. 0. 

Therefore, as a next step the relative influence of ef- 
fects due to the strength of the adsorption potential and 
to the resistance of helium atoms to be compressed was 
analyzed. It was found that if one modifies /3ho by in- 
cluding a factor which takes into account the ratio be- 
tween the well depth D and the compressibility at sat- 
uration density 1/poKy = 27.2 K 7] in such a way that 



Pmho = /9ho 



D 



1/4 



3 m tr?. 



/ D 



1/4 



PO 



1/4 



(26) 



one gets a good overall agreement as shown in Fig.[S| It is 
worthy of notice that the slope of this simple expression 
is written in terms of physical observables of bulk helium. 

Moreover, it is important to emphasize that Fig. 
collects data on substrates from H2 and Li, which are 
barely able to support monolayer helium films, up to 
Gr (graphite), which is the strongest attractor for he- 
lium atoms Equation H26|) can be considered as a 
nearly universal property since the straight line for /3mho 
matches well the data of /3pH for all the substrates in- 
cluded in the present Table HI i.e., for those plotted in 
Fig. 2 of Ref. yj provided they ly above the dashed curve 
traced there. Notice that in that figure the points deter- 
mined by D and C3 are spread over the whole drawing. 

In conclusion, one can state that a simple semi- 
phenomenological correction in the surface term of the 
SFM given by Eq. Q yields an important improvement 
of the prediction for the E ^ F phase transition. The 
corrective term proposed in the present report contains a 
decaying exponential built by following the simple model 
formulated by Cheng et al. The cutoff parameter 

can be interpreted by taking into account properties of 
the harmonic g.s. wavefunction, the well depth of the 
adsorption potential and the compressibility of bulk he- 
lium. 
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